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Abstract: [ Aim] In this study, the potential side effects of the sublethal doses of imidacloprid on the 
midgut bacteria of the worker bees of Apis mellifera ligustica was investigated under laboratory conditions 
so as to provide a theoretical basis for the biosafety assessment of pesticides. [ Methods] Newly emerged 
bees were fed on syrups with different sublethal doses of imidacloprid (0.005 mg/L, 0.015 mg/L and 
0.045 mg/L) and pure sugar syrup with 0.3% acetone, and samples were collected after 3, 6, 9, 12 
and 15 d, respectively. The population structure of bacteria in the midgut of honeybees in the four 
treatment groups was analyzed using next-generation sequencing after the genomic DNA was gained from 
intestinal bacteria. 【 Results] The results of sequence analysis showed that the midgut bacteria in A. 
mellifera ligustica mainly belong to Proteobacteria and Firmicutes. The composition of dominant midgut 
bacteria and the richness estimators in the same period were analyzed for statistical significance among the 
four treatments by performing a T-test, and no significant differences were observed ( P > 0.05). 
【 Conclusion] The results suggest that the sublethal doses of imidacloprid does not induce significant 
changes in the composition of the honeybee gut bacterial community under laboratory conditions. 
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1 INTRODUCTION 


Neonicotinoids are a group of the most broadly 
used insecticides with more than 25% share of global 
market (Van der Sluijs et al., 2013). Imidacloprid 
is a representative of neonicotinoid insecticide. Due 
to its high efficiency, low toxicity, and broad- 
spectrum, imidacloprid has become the second most 
used agrochemical in the world. 

Honey bees are the most important pollinators in 
and play an 
in biodiversity maintenance and 


natural ecosystems world-wide, 
important role 
ecosystem function. They, therefore, have been 
widely considered as an important insect in the 
biosafety assessment of insecticides. The population 
number of honeybees has declined in recent years 
(Brown and Paxton, 2009; Cameron et al., 2011). 
Many investigators suspected that extensive pesticide 
use is one of the major contributors to heavy bee 


losses worldwide ( vanEngelsdorp and Meixner, 
2010; Blacquiere et al., 2012). In May 2013, the 
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adopted a 2-year 
restriction to the use of three neonicotinoid pesticides 


European Commission ( EC ) 


(imidacloprid, clothianidin, and thiamethoxam) on 
crops that attract bees and other pollinators in an 
effort to reduce bee losses and to re-evaluate the risk 
of these pesticides on honeybees within the two years 
(Maxim and Arnold, 2013). 

Many environmental factors affect bees, but do 
not cause direct mortality. Those factors are called 
sublethal impact. Honeybees tend to be exposed to 
pesticides at sublethal doses in the real situation, so 
assessment of sublethal effects of pesticides on bees 
has been a subject with increasing discussion within 
scientists ( Laurino et al., 2011). Sublethal effect of 
pesticides on honeybees has been most widely 
studied with the focus mainly on foraging behavior 
2010 ), 


performance and memory ( Decourtye et al., 2011) , 


( Decourtye and Devillers, learning 


motor activity ( Lambin et al., 2001) , and olfactory 
capacity (Hassani et al., 2005). 

A number of findings are revealing that there is a 
between the intestinal 


strict interconnection 
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microbiota balance and the health status of the 
honeybees ( Cox-Foster et al., 2007; Fraune and 
Bosch, 2010; Engel et al., 2012). The possibility of 
using the gut bacterial community as an indicator of 
honeybee health has been discussed ( Crotti et al., 
2013), but there are no reports focusing on the 
midgut bacterial structure in honeybees exposed to 
imidacloprid. The aim of the present study is to 
investigate the impact of imidacloprid, administered 
at sublethal doses, on the midgut bacterial community 
of worker bees over time under laboratory conditions. 


2 MATERIALS AND METHODS 


2.1 Insect and sublethal doses of imidacloprid 
preparation 

Newly emerged worker bees (less than 1 h old) 
were chosen in the experiment; the bees were 
randomly sampled from the apiaries of Department of 
Bee Protection and Biological Safety, Institute of 
Apicultural Research, Chinese Academy of 
Agricultural Sciences, Beijing, China. 


Imidacloprid (99.99% pure ) 
powdered form from Sigma-Aldrich (Germany). To 


obtained in 


identify whether the concentration of imidacloprid 
influenced the intestinal bacterial community of the 
honeybees, three concentrations (0. 005, 0. 015, 
and 0.045 mg/L) of imidacloprid were used in the 
experiment. A stock solution of imidacloprid was 
prepared in acetone. Aliquots of the stock solution 
were used to make each test solution at a specific 
concentration, which were added to sugar syrup 
[50% (w/v) sucrose solution | to obtain the desired 
imidacloprid concentrations. The final concentration of 
acetone in the sucrose solutions was 0.3%. The effects 
of imidacloprid-containing solutions were compared with 
the control group that received the treatment of sucrose 
solution containing 0.3% acetone. 
2.2 Laboratory feeding of honeybees 

Newly emerged bees (less than 1 h old) were 
assigned randomly to wooden cages ( dimensions of 
10 cm x7 cm x8 cm) with mesh on two sides. For 
each treatment, a total of 80 bees were used. All the 
cages were kept in a dark incubator at 29 + 1°C and 
65% + 5% 
workers possess few bacteria ( Martinson et al., 
2012), so bees were fed with 2 000 pL of sugar 
syrup [50% (w/v) sucrose solution], 300 uL of 
water ad libitum, and 0.5 g of pollen from a gravity 


relative humidities. Newly emerged 


feeder fitted on each cage for the first 7 days 
(Ramirez-Romero et al., 2008; Dai et al., 2012). 
From the eighth day, the laboratory honeybees were 
exposed to sublethal doses of imidacloprid, and 33 
ML of the sugar solution per bee was supplied every 


day (Decourtye et al., 2003). Four treatments were 
designed to examine the effects of the sublethal doses 
of imidacloprid on the intestinal bacterial community 
of the honeybees, and for each treatment, three 
replicates were undertaken. A group of bees receiving 
the treatment of sucrose solution containing 0. 3% 
a negative control group. 
Honeybees of the other three treatments were fed on 
syrups with 0.005 mg/L, 0.015 mg/L, and 0.045 mg/ 


L of imidacloprid, respectively. Bees were collected for 


acetone is served as 


further analysis following exposure to imidacloprid for 
3, 6,9, 12 and 15 d, respectively. 
2.3 Honeybee dissection and DNA extraction 

The honeybees were dissected at 3, 6, 9, 12 
and 15 d after treatment, respectively. The bees 
were placed at — 20°C for 20 s before dissection. 
The midguts of the honeybees were isolated in sterile 
conditions, and ten samples were collected for each 
treatment, kept in 2. 0 mL eppendorf tubes and 
- 70°C until use. DNA 
extraction was performed according to the QIAamp 
DNA Stool Mini Kit protocol (QIAGEN, Germany). 
2.4 16S rRNA gene V3 region PCR and 
pyrosequencing 

Amplification of partial bacterial V3 region of 
16S rRNA gene sequences were performed as 
previously described (Nakatsu et al., 2000) with the 
primer pair V3-343f and V3-534r. PCR reactions 
were performed in a volume of 50 uL containing 22 
uL sterile water, 1.0 pL DNA template, 1.0 uL of 
each primer, and 25 uL 2 x Master Mix ( Promega, 
USA). The following touchdown PCR protocol was 
used; initial denaturation at 94°C for 5 min, followed 
by 10 cycles at 94°C for 30 s, 60°C for 45 s, and 


72°C for 90 s, in which the annealing temperature was 


immediately stored at 


reduced by 0. 5°C/cycle from the preceding cycle, 
followed by 27 cycles of 94°C for 30 s, 55°C for 45 s, 
and 72°C for 90 s. The PCR products (about 200 bp) 
were separated by electrophoresis on agarose gels 
(1.5%, w/v) and recovered with a DNA Gel 
Extraction Kit (AxyGEN, USA). The amplified PCR 
products were sequenced by the Illumina HiSeq 2500 
according to the manufacturer’ s instructions in 
Beijing Computer Center (Beijing, China). 
2.5 Sequence Analysis 

Illumina reads were binned according to index 
sequence. The data derived from the sequencing was 
processed using FastQC program ( http://www. 
bioinformatics. babraham. ac. uk/projects/fastqc/ ) . 
After initial trimming, we obtained the merged reads 
with FLASH vl. 2. 7 ( http://ccb. jhu. edu/ 
software/FLASH/) based on overlapping regions 
within paired-end reads ( Magoc and Salzberg, 
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2011), and then removed the short sequences ( < 
150 bp) with NGSQCToolkits2. 3 (Patel and Jain, 
2012). Poor-quality sequences (more than 80% of 
the bases demonstrated a quality score of <20 in the 
raw reads ) were removed using QIIME1.8 
(Caporaso et al., 2010) , which 
were formed during PCR reaction, were discarded 
with USEARCH v6.1 (Edgar, 2010). All Illumina 


sequence data of this study have been submitted to 


and the chimeras, 


the National Center for Biotechnology Information 
(NCBI) Sequence Read Archive (SRA) and given 
the accession number SRP047113. The resultant 
high-quality sequences were clustered to operational 
taxonomic units (OTUs) at the cut-off of 97% 
similarity using UCLUST vl. 2. 22. The taxon 
abundance of each sample was generated into 
phylum, class, order, family, and genera levels 
using the RDP classifier and an 80% confidence 
threshold ( Wang et al., 2007). The Rarefaction 
curves, sample coverage and richness estimators 
were calculated using QIIME toolkit, and then the 
richness estimators were analyzed for statistical 
significance by performing a T-test. 


3 RESULTS 


3.1 Basic statistics of 16S rRNA gene sequences 
from A. mellifera ligustica 
The DNA sequences of PCR products were 
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sequenced by the Illumina HiSeq 2500. A total of 
27 673 383 of V3 16S rRNA sequences reads from 
the 60 samples with an average of 461 223 sequences 
reads for each sample (the minimum of one sample 
was 356 486 and the maximum was 529 077) were 
used for this project. The average length of sequence 
reads was 131 bp. The richness estimators and 
rarefaction curves for all samples plateaued. This 
suggests that the current analysis had already 
captured most microbial diversity, although more 
phylotypes found by 
sequencing depth. 
3.2 Intestinal bacterial communities of adult 
honeybees of A. mellifera ligustica 

Following quality filtering, we generated a data 
set consisting of 15 102 475 high-quality and 
classifiable 16S rRNA gene sequences for further 


may be increasing the 


analysis. They were classified into different taxonomy 
using QIIME. The taxon abundance of each sample 
was generated into phylum, class, order, family, and 
genera levels using mainly database of RDP, aided by 
Greengene databases. At the 
Proteobacteria and Firmicutes were the most dominant 
bacterial phyla, which occurred in 
bacterial community of honeybees ( Disayathanoowat et 
al., 2012; Geng et al., 2013). At the genus level 
(Fig. 1), Kingella of the Beta-Proteobacteria was the 


phylum level, 


consistently 


most abundant genus, and the amount of Lactobacillus 


was the second. 
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Relative abundance of the dominant bacterial genera in the midgut of Apis mellifera ligustica adults 


A-E; Exposed to imidacloprid for 3, 6,9, 12, and 15 d, respectively; Ck; Pure sugar syrup with 0.03% acetone; PS; 0.005 mg/L imidacloprid 
syrup; P15; 0.015 mg/L imidacloprid syrup; P45; 0.045 mg/L imidacloprid syrup. 
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3. 3 Effects of the sublethal doses of 
imidacloprid on the midgut bacterial composition 
of A. mellifera ligustica 

The four feeding treatments were 0.005 mg/L 
imidacloprid syrup, 0. 015 mg/L 
syrup, 0. 045 mg/L 
pure sugar syrup with 0. 3% 


imidacloprid 
syrup, and 
acetone. The 


imidacloprid 


histograms of the dominant midgut bacterial genera 
were generated from adult honeybee samples 
collected 3, 6, 9, 12, and 15 d after feeding on 
syrups with different concentrations of imidacloprid 
and pure sugar syrup (Fig. 2). In order to explore 
how the sublethal doses of imidacloprid affect on the 
composition of midgut bacteria, the composition of 
dominant midgut bacterial genera were analyzed for 
statistical significance among the four treatments in 
the same period by performing a T-test, and no 
significant differences were observed ( P >0. 05). 
This result suggests that the sublethal doses of 
imidacloprid had no effect on midgut bacterial 
composition. 
3. 4 Effects of the sublethal doses of 
imidacloprid on the midgut bacterial structure of 
A. mellifera ligustica 

The resultant high-quality 


sequences were 
clustered to operational taxonomic units (OTUs) at 
the cut-off of 97% similarity using UCLUST 
vl.2.22, then were used to calculate the 
Rarefaction curves, sample coverage and richness 
estimators in QUME platform, and these histograms 
of Shannon estimator were constructed using QIIME 
toolkit (Fig. 3), providing a clear visualization of 
the relationships among the different sample groups. 
The richness estimators, including Shannon, chaol , 
and ACE, were analyzed for statistical significance 
among the four treatments in the same period by 
performing a T-test, and no significant differences 
were observed ( P >0.05). This result suggests that 
the sublethal doses of imidacloprid had no effect on 
midgut bacterial structure. 


4 DISCUSSIONS 


Effects of neonicotinoid pesticides on honeybees 
have been most widely studied, and their adverse 
effects on foraging behavior ( Decourtye and 
Devillers, 2010) , learning performance and memory 
( Decourtye et al., 2011) , motor activity ( Lambin et 
al., 2001), and olfactory capacity (Hassani et al., 
2005) have been well analyzed. Some metagenomic 
and experimental studies suggest that gut bacteria of 
honeybees play a major role in neutralization of 


dietary toxins, nutrition, and in defense against 


pathogens (Engel et al., 2012). The possibility of 
using the gut bacterial community as an indicator of 
honeybee health was discussed ( Crotti et al., 
2013). However, the effects of the neonicotinoid 
pesticides on gut bacterial communities of honeybees 
have not been investigated yet. In this study, the 
effects of the sublethal doses of imidacloprid on the 
midgut bacterial communities of A. mellifera ligustica 
were first investigated under laboratory conditions. 
But no significant differences were found in midgut 
bacterial structures and composition among adult 
honeybees fed on syrups with different sublethal 
doses of imidacloprid (0. 005, 0. 015, and 
0.045 mg/L) and pure sugar syrup with 0. 03% 
acetone as a negative control group under the 
laboratory conditions. Moreover, the samples were 
collected at different time (3, 6,9, 12, and 15d), 
and no obvious differences were observed. 

In this study, we used Illumina next-generation 
sequencing which is one of the non-culture based 
methods, to analyze the midgut bacterial community 
of the honeybees. The microbial community can be 
evaluated with high accuracy, low cost and short 
time by using the new sequencing technologies. The 
new deep sequencing methods provide a convenient 
platform for characterizing the composition of the 
animal microbiota. Furthermore, compared with 
earlier culture-based methods, using these methods 
can retrieve the most abundant members, so these 
new techniques have been widely used in analyzing 
the animal microbiota. Also the non-culture based 
methods have been used in characterizing the 
honeybee microbiota since 2003 (Jeyaprakash et al., 
2003), and recently these techniques have been 
applied in more and more experiments on honeybee 
microbiota ( Cox-Foster et al., 2007; Engel et al., 
2012; Martinson et al., 2012). 

The dominant midgut bacterial taxa of A. 
mellifera ligustica detected in this study were B-, y- 
Proteobacteria and Firmicutes, consistent with the 
previous reports on the intestinal microbiota diversity 
of A. mellifera ( Martinson et al., 2011; 
Disayathanoowat et al., 2012; Geng et al., 2013). 
This result supports 
honeybees have a co-evolved symbiotic relationship 


the previous opinion that 


with a few specific bacterial groups once again 
(Martinson et al., 2011). However, the diversity 
found was relatively low as compared with these 
studies. Furthermore, in this study, members of the 
8-Proteobacteria were found to be the most common 
group of bacteria in honeybee midgut, which was 
contrary to some studies ( Babendreier et al., 2007; 


Geng et al., 2013). Geng et al. (2013), for example, 
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Fig. 2 Relative abundance of the dominant midgut bacterial genera in the midgut of 
Apis mellifera ligustica adults after exposure to imidacloprid for different time 
A-E; Exposed to imidacloprid for 3, 6,9, 12, and 15 d, respectively. Ck; Pure sugar syrup with 0.03% acetone; PS; 0.005 mg/L imidacloprid 
syrup; P15; 0.015 mg/L imidacloprid syrup; P45; 0.045 mg/L imidacloprid syrup. 
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Fig. 3 Shannon index of bacterial diversity in the midgut of Apis mellifera ligustica adults after exposure 
to imidacloprid for different time calculated using QUME toolkit 
A-E; Exposed to imidacloprid for 3, 6,9, 12, and 15 d, respectively. Ck; Pure sugar syrup with 0.03% acetone; P5: 0.005 mg/L imidacloprid 
syrup; P15; 0.015 mg/L imidacloprid syrup; P45; 0.045 mg/L imidacloprid syrup. 


reported that y-Proteobacteria was the most common 
group of bacteria in honeybee midgut. Lactobacillus is 
a major diversity genus in this study, which is in 
accordance with the previous studies ( Vasquez et al., 


2012; Disayathanoowat et al., 2012). 


evidence shows that Lactobacillus plays a significant 


Increasing 


role in maintaining and improving honeybee health 
(Evans and Lopez, 2004). 

Our study provides an insight into the potential 
risk of imidacloprid on the bacterial diversity in the 
midgut of A. mellifera ligustica. It concludes that 
sublethal doses of imidacloprid do not impact the 
bacterial diversity in the midgut of honeybees (A. 
mellifera ligustica ) laboratory 
Therefore , total bacterial community structures may not 


under conditions. 


be a sensitive indicator for providing evidence for the 
impact of pesticide on honeybees at sublethal levels. 
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亚 致 死 剂量 吡虫啉 对 意大利 蜜蜂 
中 肠 细 画 群 洛 的 影响 


(中 国 农业 科学 院 蜜蜂 研究 所 ,农业 部 授粉 昆虫 生物 学 
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摘要 :【 目 的] 在 实验 室 
以 期 为 杀 虫 剂 的 生物 安 
的 吡虫啉 蔗糖 溶液 以 及 
12 和 15 d 时 分 别 取样 。 























含 丙酮 (0.039%6 ) h 








条 件 下 ,研究 亚 致死 剂量 吡虫啉 对 意大利 蜜蜂 Apis mellifera ligustica 中 肠 细菌 群落 的 影响 ， 























全 性 评价 体系 提供 理论 依据 。[ 方 法 ] 分 别 用 3 种 不 同 浓度 (0. 005, 0.015 和 0.045 mg/L) 

















提取 中 肠 细菌 基 

















9 糖水 (空白 对 照 ) 在 实验 室内 饲 喂 刚 出 房 的 意大利 蜜蜂 ,在 饲 喂 3,6, 9， 



































因 组 DNA ,采用 下 一 代 高 通 量 测序 技术 对 4 个 处 理 组 蜜蜂 中 肠 细菌 群落 





























进行 检测 分 析 。[ 结果 细菌 序列 分 析 表 明 , 意 大 利 蜜蜂 中 肠 细菌 主要 属于 变形 菌 门 ( Proteobacteria) 和 厚 壁 菌 站 


( Firmicutes) ,同一 时 间 内 各 处 理 组 间 主 要 细菌 数量 差异 不 显著 (P >0.05) ,同时 细菌 群落 多 样 性 指数 差异 也 不 显 


著 (P>0.05)。[【 结 论 ] 在 实验 室 条 件 下 亚 致 死 剂 量 吡虫啉 对 意大利 蜜蜂 中 肠 细菌 群落 结构 影响 不 显著 。 
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